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ABSTRACT 


Y 

A  dynamic  structural  model  of  a  submerged  ring  is 
developed  using  trigonometric  series.  It  is  constructed 
for  use  in  conjunction  with  a  finite  element  fluid  model 
to  examine  the  effects  of  cavitation  on  underwater  shock 
loading  of  a  structure.  The  governing  equations  and  the 
time  integration  algorithm  used  in  the  model  are  described. 
Results  predicted  by  the  model  are  compared  to  known 
results.  The  program  listing  is  given. 
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I .  INTRODUCTION 


The  program  evolved  during  this  study  models  a  submerged 
circular  ring  using  trigonometric  series.  It  was  developed 
for  use  in  conjunction  with  a  finite  element  fluid  model 
based  on  a  displacement  potential  formulation.  The  purpose 
of  the  combined  models  is  to  predict  the  effects  of 
cavitation  on  underwater  shock  loading  of  the  structure. 

More  information  on  the  fluid  model  can  be  found  in 
References  1  and  2.  The  purpose  of  this  paper  is  to 
describe  the  development  of  the  structural  model  and  to 
present  some  of  the  results  obtained  from  its  use  in 
combination  with  the  fluid  model.  A  listing  of  the  FORTRAN  IV 
program  implementing  the  structural  model  is  given  in  the 
Appendix. 
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II.  GOVERNING  EQUATIONS 


A.  DIFFERENTIAL  EQUATIONS  FOR  THE  DEFLECTION 
OF  A  THIN  CIRCULAR  RING 

In  Figure  1  the  solid  line  represents  the  ring  after 
deformation,  and  the  dotted  line  the  ring  before  deformation. 
For  small  deflections  the  curvature  of  the  element  m^n^  can 
be  taken  as 


_1_  =  d9  +  Ad  9 
R^  ds  +  Ads 


where  R,  w,  9  and  s  are  defined  as  shown  in  Figure  1  and 
w  is  taken  positive  inward.  Making  use  of  the  relations 


Ad9  =  d9^  -  d9 


dw 

ds 


ds 


Ads  =  ds^  -  ds 


(r-w)d9  -  rd9  =  -wd9 


-w 


ds 

R 


and  substituting  into  the  equation  above  yields 


_L_ 

R, 


d9  +  ds 


ds 

dsd-5) 


d9  (1  +  |) 
ds"(1.w)(1  +  w) 


,2 

d  w  .  w. 

TT  ds(1  +  r 

ds _ 

ds(l-g)  (1  + J) 


Neglecting  the  higher  order  terms,  the  above  expression 
becomes 
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For  a  ring  where  the  thickness  is  small  compared  to  the 
radius  and  elastic  behavior  is  assumed,  it  can  be  shown 
that  the  approximate  relationship  between  deflection  and 
loading  is  [Reference  7] 

I 


where  M  is  the  bending  moment  about  the  centroidal  axis 
and  D  is  the  flexural  rigidity  of  the  ring.  A  positive 
bending  moment  produces  compression  in  the  outside  fibers 
of  the  ring.  Combining  equations  1  and  2  yields  the 
differential  equations  for  the  deflection  of  the  ring  given 
below. 


B.  GOVERNING  EQUATIONS  OF  MOTION 


The  governing  equations  of  motion  used  in  the  program 
were  arrived  at  by  the  application  of  Hamilton's  principle. 

In  order  to  apply  Hamilton's  principle,  it  is  first  necessary 
to  derive  expressions  for  the  strain  energy  and  kinetic 
energy  of  the  ring  as  well  as  the  work  done  by  the  external 
loads.  In  these  derivations  the  ring  was  taken  to  be 
homogeneous,  elastic,  and  of  unit  width.  The  pressures  on 
the  ring  and  its  deflection  were  represented  by  the 
pressures  at  and  the  deflection  of  a  set  of  nodal  points 
equally  spaced  along  the  circumference  of  the  ring. 

The  shock  front  is  assumed  to  approach  the  ring  normal 
to  the  0=0  plane  as  shown  in  Figure  2,  where  9  is  taken 
to  be  positive  counterclockwise. 

K - 


Figure  2.  Ring,  Shock  Front  Orientation 
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Since  the  ring  is  symmetric  about  the  9=0  plane, 
results  are  given  for  9  from  0°  to  180° .  The  response  of 
the  ring  is  represented  by  the  response  of  a  set  of  N+l 
nodal  points  equally  spaced  along  the  ring.  The  radial 
displacement  is  approximated  by  a  trigonometric  cosine 
series 


N 

w  =  l  a  cos  n9  (4) 

n=0 


The  tangential  deflection  is  represented  by  the  trigonometric 
sine  series 


N 

v  =  l  b  sin  n9  ,  (5) 

n=l  n 

where  v  is  taken  to  be  positive  in  the  negative  9  direction. 
Similarly  the  normal  pressure  applied  to  the  ring  is 
represented  by 


N 

p  =  £  c  cos  n9  .  (6) 

n=0  n 

1.  Strain  Energy 

The  strain  energy  is  comprised  of  two  components 
[Ref.  4] .  One  component  is  the  strain  energy  due  to  bending 
and  the  other  to  strain  in  the  9  direction  (eQ).  The  strain 
energy  stored  in  the  ring  due  to  bending  is  calculated  as 
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the  work  done  by  the  internal  moment  acting  over  the 
entire  ring,  and  for  a  ring  of  unit  width  is 


UB  =  I  /  M  R<de 
0 

where  <  is  the  curvature  of  the  ring.  Combining  this 
equation  with  equation  3  and  using  the  relation  <  =  M/D 
yields 


u. 


B 


"  J  / 


2  7T 


1 

? 


~2 
-  3  w 

2 

39^ 


+  w) 


de 


(7) 


The  strain  energy  due  to  the  strain  in  the  9  direction  is 
obtained  as  the  distance  traveled  by  the  average  normal 
force  in  the  9  direction. 

The  circumferential  strain  is  composed  of  two 
components  the  first  of  which  is  the  result  of  the  radial 
displacement  of  the  ring  (  )  and  the  second  results  from 

tangential  displacement  and  is  -|^)  .  The  strain  energy 
from  the  normal  force  (N)  is 


0, 


-  7  / 


2  TT 


neqr  de 


"  T  / 


2n 


ASqEEqR  d0 


or 


(8) 


Ur 


-  J  / 


2  it 


RAEe, 


de  = 


1  /'  AZ 

2  J 


j  <U +  ">2  Me 
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where  A  is  the  cross  sectional  area  of  the  ring  and 

E  is  the  modulus  of  elasticity  of  the  ring  material. 


1 


Combining  equations  7  and  8,  the  expression  for  the 


total  strain  energy  is 


I  /2"4  +  I  /*'  ¥‘S*«>2« 

o  R  39  o 

o  o  (9) 


Substituting  into  this  equation  the  series  expressions  for 
radial  and  tangential  displacements,  and  assuming  constant 
geometrical  and  material  properties  over  the  ring  yields 
after  integration 


N  M  * 

„  tD  ,  ?  2.  2  ,r  AEtt  ,  v  ,  .  .  ,2, 

UT  =  —  1  I  an  (n  "  1)  +  ~2R  C  I  (nbn  +  an}  1 

1  2R  J  n=0  n  n=0  n  n 


An  expression  for  the  change  in  strain  energy  SU 
corresponding  to  small  changes  in  tangential  and  radial 
displacements  (6i^and5an)  can  now  be  found 

N  N 

<5U  =  £2  [  I  (n2  -  1)  2a  6a  ]  +  ^  [  £  (n2b  +  na  )  6b 

3  n=0  nn  R  nilQ  n  nn 

+  (nb  +  a  ) 6a  ) ]  (10) 

n  n  n 


In  this  equation  and  several  that  follow  which  contain 
a  summation  expression  starting  with  n  =  0,  a  factor  of  2 
associated  with  the  n  =  0  term  is  omitted  for  the  sake  of 
brevity.  The  term  is  accounted  for  in  the  solution  process 
used  with  the  model. 


18 


2.  Work 


The  work  associated  with  a  small  change  in  radial 
displacement  is  given  by 


2tt 

5W  =  /  p  5w  R  d0 

0 


Substituting  the  relations 


N 

p  =  V  c  cos  m0 
r  L  n  m 

m=0 

and 

N 

5w  =  F  5a  cos  n0 
_  «  n 


into  this  expression  yields,  after  integration 


N 

6W  =  7  c  5a  itR 

_ ^  ^ 


3.  Kinetic  Energy 

The  expression  for  the  kinetic  energy  of  the  ring 
can  be  arrived  at  by  considering  an  infinitesimal  section 
of  the  ring  that  has  been  set  in  motion.  The  kinetic  energy 
of  the  element  can  be  represented  as  the  sum  of  a  contribu¬ 
tion  due  to  translation  of  the  mass  center  and  a  contribution 
due  to  rotation.  This  leads  to 
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dT  =  i  p  (ARd9 )  (w^  +  v*)  +  J  pIRd9(^-) 


where  I  is  the  area  moment  of  inertia  about  the  centroidal 


axis . 


Using  the  series  relations 


N  N 

w  =  T  a  cos  n9  w  =  T  a  cos  n6 

n  L  _  n 

n=0  n=0 


N 

r» 

v  =  l 

n=  1 


yields 


1  N  .  2  N  •  2 

dT  *  ■=■  p(ARd9)[(  l  a  cos  n0)  +  (  l  b  sin  n6 )  ] 

1  n=0  n=l 


+  j  pIRd9 (-  ^  l  an  n  sin  n9)z 
n=l 


Integration  over  the  ring  yields 


T  =  ^  l  [(AR  +  n2  |)an2  +  ARt>n2] 

n=0 


The  resultant  change  in  kinetic  energy  due  to  a 
small  change  in  velocity  components  can  then  be  found  as 


6T  =  TTP  l  [  (AR  +  n2  §)an6an  +  ARb  6bn]  (12) 

n=0 


N 

b  sin  n9  v  =  7  fc>  sin  n9 

n  n=l  n 
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4 .  Application  of  Hamilton's  Principle  To  Find 
Coupled  Equations  of  Motion 

From  Hamilton's  principle  [Ref.  5]  it  is  known 


that 


/  (<5T  -  5U  +  5W)  dt  =  0 


(13) 


Combining  equations  10,  11,  12  and  13  yields,  after 
carrying  out  an  integration  by  parts  on  the  first  term, 


N 


up  (AR  +  n2  i)  [a  6a 
n  R  n  n 

n=0 


2  .. 

/  an  dt] 


N 

+  l  7rpAR[b  5b 
n  n  n 

n=0 


2  .. 

/  b  6b  dt] 
'  n  n 


C2 ,tD  ,  ?  ,_2  ,,2 


-  n?'  I  (n  *1)  a„  «a„i  * 

R  n=0 


AEir 

R 


N 


[  l  (n  b  +  na  >  5b 
n=0 


N 

+  (nb  +  a  ) 6a  ]  -  7  c  6 a^  ttR}  dt  =  0 

n  n  n  n  n 

n=u 


Since  the  6a  and  5b  can  be  chosen  arbitrarily,  this 
n  n 


becomes 


...  ,r.  2  2,”  ,  ,c.2r.  ,r,  2,  2  .,2.  ,c.  2  n 

(1  +  (R>  n  an  +  (R>  [1  +  R  (n  _1)  an  +  (R}  pA 


il  .  ,c,  2  ,  ,c.  2  2,  n 

b  +  (=)  na  +  (p)  n  b  =  0 

n  R  n  R  n 


(14) 


where  c  =  (E/p)1^2  is  the  speed  of  sound  in  the  ring  and 
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where  r  is  the  centroidal  radius  of  gyration  of  the  ring 
section  and  the  relation  I  =  Ar2  has  been  used. 

5 .  Separated  Equations  of  Motion 

It  is  common  practice  to  separate  ring  deformations 
into  flexural  modes  and  extensional  modes.  The  initial 
version  of  this  program  utilized  such  a  resolution.  The 
flexural  mode  is  characterized  by  the  requirement  that 

3v 

£e  »  -  <w  +  Je)/R  =  0  • 

This  implies  that  the  Fourier  coefficients  satisfy  the 
relation 


n(F) 


+  nb 

n 


(F) 


0 


(15) 


The  extensional  mode  is  defined  to  be  geometrically  orthogonal 
to  the  flexural  mode  so  that 


na  -  b 
n(E)  n(E) 


=  0 


(16) 


Substituting  equation  15  into  equations  14  yields  the  result 


[pA(  (—  Tj-— )  +  zn2)]a  +  z  (n2  -  1)  2]  a 

n  n(F)  R2  n(F) 


=  c 


n 


5  =  -a  /n 

(F)  n(P) 


(17) 


where  z  =  (r/R) 
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If  equation  16  is  substituted  into  equations  14,  the  result 
is 


tpA  (  (n2  +  1)  +  zn2 )  ]  a  +  ((n2  +  l)2  +  z(n2-l)2]an 

n(E)  R2  n(E) 


=  c 


n 


(18) 


n 


=  na 


(E) 


n 


(E) 


Initial  solutions  were  carried  out  using  equations 
17  and  18  and  then  combining  the  results  by 


a 


n 


a  +  a 

n(F)  n(E) 


b 


n 


b  +  b 
n(F)  n(E) 


This  procedure  gave  results  which  satisfactorily  predicted 
the  ring  bending  moments  but  gave  poor  accuracy  for  axial 
force  determinations. 

Re-examination  of  the  foregoing  solution  process 
revealed  that  the  mode  shapes  defined  by  equations  15  and 
16  are  not  orthogonal  with  respect  to  the  mass  or  the 
stiffness  matrices  of  the  system.  It  was  accordingly 
decided  to  return  to  equations  14  and  solve  them 
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simultaneously  in  the  time  integration  algorithm.  This 
change  led  to  accurate  axial  force  determinations.  The 
separated  equations  of  motion,  17  and  18,  are  used  in  the 
program  as  described  in  Section  III,  b. 
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III.  TIME  INTEGRATION  METHOD 


A.  CENTRAL  DIFFERENCE  ALGORITHM 

Time  integration  is  accomplished  using  a  central 

difference  algorithm.  If  a  ,  b  ,  and  c  are  known  at 

n  n  n 

time  t^,  then  a^^  (the  value  of  an  at  t^)  may  be 
evaluated  from  the  first  of  equations  14.  Letting  h 

represent  the  time  step  and  a  ^  ^  represent  an  at  time 

( i )  * 

t  -  h/2,  the  next  value  of  an  is  calculated  from 


• (i+4) 

a 

n 


+  ha 


(i) 

n 


(19) 


Using  this  result,  the  next  an  (at  time  t^+1^  =  t^  +  h) 
is  calculated  from 


(i+l) 

u 

n 


+  h  a 


(i+*s) 

n 


(20) 


The  value  b^  is  similarly  found  from  the  second  of 
n 

equations  14.  A  pair  of  equations  paralleling  19  and  20 

are  used  to  calculate  b  and  b^.  When  these  steps 

n  n 

have  been  completed  for  each  n,  the  value  of  radial 
displacement  w  is  found  at  each  structural  node  from 
equation  4.  These  displacements  are  passed  to  the  fluid 
program  where  they  furnish  required  interface  boundary 
conditions  to  allow  an  advance  to  time  t1''*'4"'^  .  Values  of 
fluid  pressure  at  the  interface  nodes  are  returned  by  the 
fluid  program. 
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The  new  nodal  pressures  are  used  to  calculate  cn ' s  at 

t  ( i-*"! )  .  it  is  then  again  possible  to  advance  an  and  an 

using  equations  14,  19,  and  20  and  to  perform  the  parallel 

calculations  for  b  and  b  . 

n  n 

The  solution  process  is  started  with  the  ring  at  rest 

under  uniform  hydrostatic  pressure.  Under  the  loading  a^ 

is  the  only  nonzero  Fourier  coefficient.  Because  a  ^  =  0 

n 

.  t-k) 

is  known  (rather  than  a^  ,  a  fictitious  starting  value  of 
.  (-k) 

a.  'is  first  calculated  from 
n 


(-H)  _  h  "  ( 0 

n  "  "  2  a 


A  similar  starting  procedure  is  used  for  bn- 

B.  SELECTION  OF  THE  NUMBER  OF 
VIBRATIONAL  MODES  USED 

The  number  of  vibrational  modes  that  can  be  modeled 
accurately  is  limited  by  the  numerical  integration  algorithm. 
Specifically,  the  algorithm  becomes  unstable  for  time  steps 
in  excess  of  about  0.3  of  the  period  of  the  structural  mode. 
The  accuracy  of  the  algorithm  deteriorates  even  before  the 
stability  limit  is  reached.  For  this  reason,  a  criterion 
was  established  for  limiting  the  number  of  modes  used, 
based  on  the  time  steps  selected. 

In  the  program  the  separated  equations  of  motion  for 
the  extensional  and  flexural  modes  were  used  to  find  the 
frequencies  needed  in  applying  the  criterion.  From 
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equations  17  and  18  the  natural  frequencies  of  the 
extensional  and  flexural  modes  were  found  to  be 


2  .c*2. i  2. 

u>  =  (p)  (1  +  n  ) 

n(E)  R 


n  =  0 , 1 , 2 , 3 ,  .  .  . 


m  2  -  /C,  2  ,r.  2  n2  (n2  -  l)2  n  _  ,  ,  , 

(i)  —  V'q  )  V  '5')  -5  n  —  2,0, 4  /  . 

n(F)  R  R  n2  +  1 


The  criterion  used  was  that  the  period  for  the  highest 
mode  retained  be  at  least  five  times  the  time  step  used. 

Since  the  extensional  modes  have  higher  frequencies  for  a 
given  n,  this  criterion  comes  into  effect  first  for  the 
extensional  modes.  Once  the  point  is  reached  (where 
n  =  nE(max)),  the  program  switches  from  the  coupled  equations 
of  motion  14 ,  to  the  separated  equation  17 .  When  this 
occurs  the  higher  extensional  mode  coefficients  (n  >  n,,(max)) 

£j 

are  not  needed,  and  only  the  flexural  mode  coefficients 

(for  n  >  n_ (max) )  are  used.  If  the  time  step  is  large 

enough  that  the  criterion  is  also  met  by  the  flexural 

modes  (n  =  n- (max) )  then  both  the  a  ’  s  and  b  's  are  omitted 
f  n  n 

for  all  modes  where  n  is  greater  than  n_  (max) . 

r 
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IV.  RESULTS 


To  check  the  accuracy  of  the  results  obtained  from  the 
program,  two  kinds  of  comparisons  were  made.  First,  the 
solution  of  equations  14  arrived  at  by  the  program  were 
compared  to  an  analytic  solution  of  equations  14.  Second, 
to  verify  that  equations  14  correctly  predict  ring  behavior, 
calculated  dynamic  response  of  the  ring  to  two  particular 
loadings  was  compared  to  known  static  results  [Reference  6] 
for  the  same  loadings. 

A.  COMPARISON  OF  MODEL  RESULTS  TO  EXACT  SOLUTION 

The  check  on  solution  accuracy  was  based  on  the  exact 
solution  for  a  ring  initially  at  rest  and  suddenly  loaded 
by  a  steady  pressure 

P  =  [1.625  -  .25  cos  0  -  .375  cos  29]  MPa 

The  exact  solution,  using  the  parameters  specified  in 
Table  I,  is 

aQ  =  4. 05844xlC)”3  [1  -  cos  (1011. 303t)  ] 

a1  =  -159 . 483t2  -  1 . 5617xl0*4 [1  -  cos  (1429. t) ] 

a2  =  -.10412911  -  cos  (85.616) t] 

-  3.7355xl0"5 [1  -  cos  (2259. 7t) ]  (22) 
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Table  I  gives  results  for  radial  displacements  at 
t  =  1  msec  for  9=0°,  90°  and  180°.  Computer  results  are 
given  for  four  different  time  steps.  Even  the  coarsest 
time  step  gives  results  within  1%  of  the  exact  values 
obtained  from  equations  22. 

B.  COMPARISON  OF  MODEL  RESULTS  TO  STATIC  BEHAVIOR 

To  compare  the  dynamic  response  predicted  by  the  program 
to  the  known  static  response  of  the  ring  for  the  same 
loading,  two  simple  cases  of  loading  were  assumed.  The 
first  loading  case  examined  was  uniform  pressure  surrounding 
the  ring;  the  second  case  was  two  equal  and  opposite  concen¬ 
trated  loads  acting  180°  apart  as  shown  in  Figure  3.  The 
second  loading  case  could  not  be  represented  exactly  by 
the  finite  trigonometric  series  loading  representation  used 
by  the  program.  It  was  approximated  by  nonzero  values  of 
pressure  only  at  9  =  0°  and  9  =  180°. 


B 


Figure  3.  Concentrated  Loading  on  Ring 
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TABLE  I 


Time 

Step 

sec 

AR 

NP1 

m  x  10  7 

AR 

NP2 

m  x  10  7 

E 

( 

AR 

predicted 
by  program 

2 . 5xl0~4 

11619. 

23687. 

17556. 

lxlO-5 

11679 

23472. 

17553 

lxlO-6 

11682 

23469. 

17555. 

lxlO'7 

11683. 

23469. 

17555. 

AR  frcm 
equations  17 

11683. 

23469. 

17555. 

.  2 
Ring  Parameters:  R  =  5  m,  A  =  .05  m  ,  E  =  200  GPa 

p  =  7830  kg/rn^,  r  =  .158m 


Comparison  of  program  solution  with 
analytic  solution 


TABLE  II 


AR  (max) 

mm 

Axial 

Force 

(Max) 

MN 

Program  results 

4.98 

-9.8 

TWice  Known 

5.00 

-10.0 

Static  Results 

2 

Ring  Parameters:  R  =  5  m,  A  =  .05  m  ,  E  =  200  GPa 
p  =  7830  kg/hi^ '  r  _  m 


Comparison  of  program  dynamic  solution 
with  twice  known  static  solution  for 
uniform  external  pressure  on  ring 
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Comparison  between  the  dynamic  and  static  responses 
of  the  ring  were  made  for  the  radial  displacements,  bending 
moments,  and  forces  at  particular  nodal  points  on  the  ring. 
For  the  uniform  pressure  case  any  nodal  point  on  the  ring 
can  be  used  since  the  displacements  are  uniform  over  the 
ring.  For  the  concentrated  loading  case,  the  nodal  points 
chosen  were  at  the  points  of  application  of  the  forces  and 
midway  between  the  application  points. 

For  the  dynamic  response  of  the  ring  predicted  by  the 
program,  the  ring  is  taken  to  be  initially  at  rest  and 
undeformed.  The  anticipated  behavior  of  the  ring  after 
the  application  of  the  loading  is  for  the  ring  to  deflect 
through  the  radial  displacements  found  for  the  static 
loading  and  to  continue  to  deflect  until  a  radial  displace¬ 
ment  with  a  magnitude  of  approximately  twice  the  static 
displacement  is  reached.  The  maximum  deflections  were 
expected  to  occur  after  a  time  approximately  equal  to 
one  half  the  period  of  the  vibrational  mode  that  dominates 
the  deflection  of  the  ring.  For  the  uniform  loading  case 
the  dominant  vibrational  mode  is  the  fundamental  extensional 
mode.  For  the  concentrated  loading  case,  it  is  the  second 
flexural  mode.  The  periods  for  these  vibrational  modes  can 
be  found  using  equations  21.  Tables  II  and  III  give  the 
comparisons  between  the  dynamic  response  from  the  program 
and  twice  the  known  static  response  for  the  two  loading 
cases  considered.  As  is  seen,  relatively  good  correlations 
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TABLE  III 


Slumber  of 
Stodal  Points 
Jsed 

ARA(max) 

ARg(max) 

Axial 
Force 
at  B 

Bending 
Moment 
at  B 

m 

m 

MN 

MN.m 

Model 

Results 

.2965 

-.2839 

-3.760 

-7.910 

Static 
Analysis 
Results 
Times  2 

.2634 

-.2421 

-3.536 

-6.434 

7 

Model 

Results 

.2013 

-.1857 

-2.353 

-4.737 

Static 
Analysis 
Results 
Times  2 

.1862 

-.1713 

-2.500 

-4.543 

q 

Model 

Results 

.1510 

-.1393 

-1.743 

-3.544 

Static 
Analysis 
Results 
Times  2 

.1425 

-.1311 

-1.913 

-3.477 

2 

Ring  Parameters:  R  =  5  m,  A  =  .05  m  ,  E  =  200  GPa 
0  =  7830  kq/m.3,  r  =  .158  m 


Comparison  of  dynamic  solution  with  twice  known 
static  solution  for  loading  of  Fig.  3 
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between  the  dynamic  responses  and  twice  the  static  responses 
were  obtained.  The  maximum  values  shown  in  Table  III  occur 
at  approximately  35.3  milliseconds  after  time  zero,  which 
compares  favorably  to  one  half  of  the  period  of  the  second 
flexural  mode  (36.6  milliseconds).  The  maximum  values  in 
Table  II  occur  at  3  milliseconds  which  is  approximately 
one  half  the  period  of  the  fundamental  extensional  mode 
(3.11  milliseconds). 

C .  OUTPUT  EXAMPLES 

A  sample  of  the  output  produced  by  the  structural 
program  is  shown  in  Figures  4  and  5.  The  print  code  at  the 
top  of  Figure  4  is  an  input  code  used  to  select  the  informa¬ 
tion  desired  in  the  print  out.  Following  the  print  code 
are  two  sets  of  input  parameters  to  be  used  in  the  structural 
model.  These  are  followed  by  the  mass  and  stiffness 
coefficients  calculated  for  each  of  the  vibrational  modes. 
This  information  remains  constant  for  any  given  run  and 
is  therefore  only  printed  once.  The  rest  of  the  information 
varies  with  time  and  can  be  printed  as  often  as  desired. 

The  pressures  at  nodal  points  generated  by  the  fluid  model 
are  printed  out;  these  pressures  are  for  one  time  step 
before  the  time  printed  below  the  pressures.  Following  the 
time  come  the  Fourier  coefficients  for  the  radial  displace¬ 
ments,  tangential  displacements,  and  pressures.  In 
Figure  5  at  the  top,  the  bending  moment  and  the  axial  force 
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at  each  nodal  point  are  printed  for  the  time  given  in 
Figure  4.  Tabulation  of  the  radial  displacement  follows. 

Figure  6  is  an  example  of  the  graphical  output  obtained 
from  the  fluid  model  when  the  fluid  and  structural  models 
are  run  jointly.  It  is  a  time  sequence  of  eight  plots  of 
fluid  nodal  pressures  over  the  domain.  The  plots  are  shown 
at  8  ms  intervals.  The  left  hand  edge  of  each  plot  repre¬ 
sents  the  plane  of  symmetry  and  the  top  and  bottom  rows 
represent,  respectively,  the  entry  and  exit  faces  for  the 
shock.  The  x's  on  the  lower  left  side  of  the  plots 
represent  dummy  nodes  inside  the  structure.  The  pressure 
ranges  for  the  mapping  characters  are  shown  in  Figure  7. 

In  Figure  6  the  development  of  the  cavity  can  be  followed. 

It  develops  in  the  second,  third,  and  fourth  frames, 
collapses  in  the  fifth,  and  has  vanished  by  the  sixth. 
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Fig#  4.  Output  Example 
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Fig.  5.  Output  Example  (cent.) 
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Pressure  Range  (m  pa) 
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X 

blank 

1 

3 

6 

8 

H 

A 

B 

C 

D 

E 

F 


Structure 
Cavitated 
0  -  .04 
.04  -  .09 
.09  -  .14 
.14  -  .18 
.18  -  .43 
.43  -  .77 
.77  -  1.10 
1.10  -  1.43 
1.43  -  1.77 
1.77  -  2.10 
over  2.10 


Figure  7.  Pressure  Code 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  comparisons  made  between  the  program  results  where 
convergence  of  the  solution  has  been  obtained  to  results 
obtained  by  other  methods,  show  that  the  model  satisfactorily 
predicts  the  dynamic  response  of  a  ring  to  external  pressures. 
The  ring  model  cannot  predict  exactly  the  behavior  of  a  real 
structure  such  as  a  submarine  hull.  However,  by  a  suitable 
choice  of  ring  parameters  the  response  of  the  ring  structure 
will  provide  a  good  first  approximation  to  the  response  of 
a  real  structure. 

A  useful  improvement  to  the  program  would  be  to  include 
a  section  which  would  automatically  pick  out  maximum  values 
of  bending  moments  and  axial  forces  in  the  ring  as  well  as 
when  and  where  they  occur.  This  would  allow  any  desired 
time  intervals  between  data  printout  to  be  used  without 
missing  peak  responses. 
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APPENDIX 
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